function Y = vec(X)
Y = X(:) ;
end
